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Abstract 

We study the statistical properties of geodesies, i.e. paths of minimal length, 
' (-('■ in large random planar quadrangulations. We extend Schaeffer's well-labeled 

Qh. tree bijection to the case of quadrangulations with a marked geodesic, leading 

r-| ! to the notion of "spine trees", amenable to a direct enumeration. We ob- 

^ \ tain the generating functions for quadrangulations with a marked geodesic of 

fixed length, as well as with a set of "confluent geodesies" , i.e. a collection of 
non-intersecting minimal paths connecting two given points. In the limit of 
^ . quadrangulations with a large area n, we find in particular an average number 

^ ! 3 X 2* of geodesies between two fixed points at distance i ^ 1 from each other. 

\^ I We show that, for generic endpoints, two confluent geodesies remain close to 
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each other and have an extensive number of contacts. This property fails for 
a few "exceptional" endpoints which can be linked by truly distinct geodesies. 
Results are presented both in the case of flnite length i and in the scaling limit 



[ i oc n^/^. In particular, we give the scaling distribution of the exceptional 



points. 



X. 

^ , 1. Introduction 



The study of random maps is a fundamental issue both in mathematics, where it 
raises many interesting combinatorial and probabilistic problems, and in physics, where 
maps serve as discretizations for fluctuating surfaces in various domains. Understanding 
the statistical properties of ensembles of random maps is relevant in particular in the 
context of two-dimensional quantum gravity 0] or for the description of fluid membrane 
statistics. So far, most results deal with global properties of random maps, with much 
effort made on the enumeration of families of maps. Various methods were developed in 
this context, such as Tutte's original approach through recursive decomposition or 
the more technical framework of random matrix integrals . Correlations functions of 
various observables were obtained for ensembles of maps with possibly extra statistical 
degrees of freedom such as spins or particles, but mostly in a global way consisting on 
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averaging over the positions of the points where these observables are measured (see 
0]). Much less is known on more refined local properties of maps, such as the actual 
dependence of correlations on the distance between the points where the observables 
are measured In this respect, a particularly promising approach relies on a recent 

bijective enumeration technique which consists in a coding of maps by decorated trees. 
This approach was initiated by Schaeffer who found a one-to-one mapping between, on 
the one hand, planar quadrangulations, i.e. maps with tetravalent faces only, and on 
the other hand so-called well-labeled trees, i.e. plane trees with suitably labeled vertices 
|]^,|[. His construction was then extended to various families of maps, such as maps 
with prescribed faces valences , Eulerian maps and even maps with particles or spins 
||T0[] . The main advantage of the approach is that it keeps track of the graph distance of 
vertices from a given arbitrary origin vertex, thus giving access to statistical properties 
of maps controlled by this distance. A first application concerned the so-called two- 
point function, i.e. the average number of pairs of vertices at a fixed distance from each 
other among quadrangulations of fixed area. In particular, for large maps, this two- 
point function was shown to converge to a universal scaling function characterizing the 
density of vertices at a fixed distance from an arbitrary origin in the continuous scaling 
limit of large maps []5| JTl]] . 

In this paper, we address the question of the statistics of geodesies, i.e. paths which 
join two given vertices in a random quadrangulation, and are of minimal length among 
all paths linking these vertices. Here again, we rely on an extension of Schaeffer's 
bijection using well-labeled trees to code for quadrangulations with a marked geodesic 
path. In this setting, we have access to the statistical properties of geodesies with a fixed 
given length, as well as of confluent geodesies, i.e. collections of non-intersecting minimal 
paths connecting the same two vertices. Many results can be explicited in the limit of 
large quadrangulations, i.e. when the area of the quadrangulation tends to infinity. 
Two different regimes may be studied according to whether the length of the geodesies 
remains finite or scales as the appropriate power of the diverging area. This latter 
scaling regime is particularly important to extract universal properties characterizing 
the continuous limit of large quadrangulations, a probabilistic object whose construction 
is still under investigation [|T2l - |T4| . 

The paper is organized as follows: Section 2 is devoted to the combinatorics of 
quadrangulations with marked geodesies. We first recall in Section 2.1 Schaeffer's con- 
struction which associates to each quadrangulation a well-labeled tree. We then show 
in Section 2.2 how to extend this bijection to treat the case of quadrangulations with 
a marked geodesic path. This leads us to the notion of spine tree, i.e. a well-labeled 
tree with a distinguished linear subtree (the spine) inherited from the marked path. 
Spine trees are shown to be in bijection with quadrangulations equipped with a bound- 
ary made of two geodesies. Under the condition that these two boundary geodesies 
have no intermediate contact, the configurations can be continuously transformed into 
the desired configurations of quadrangulations with a single marked geodesic path. In 
Section 2.3, we derive explicit formulae for the generating function of quadrangulations 
bounded by two geodesic paths and for that of quadrangulations with a single marked 
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geodesic of given length. The connection between the two generating functions is ob- 
tained via a simple transformation consisting in taking an irreducible part. We finally 
show in Section 2.4 how to use our basic generating functions to enumerate quadran- 
gulations with marked confiuent geodesies. We distinguish two general cases: weakly 
avoiding confiuent geodesies made of paths that can touch each other with no intersec- 
tions or strongly avoiding confluent geodesies which have no contact vertices except for 
their common first and last vertices. Section 3 deals with the statistical properties of 
geodesies of finite length i in large quadr angulations, i.e. in an ensemble of quadran- 
gulations with fixed area n in the limit n — > oo, keeping i finite. We first explain in 
Section 3.1 how our enumerative results can be translated into averages in the ensemble 
of pointed quadrangulations of fixed area n (with inverse symmetry weights) or the 
ensemble of quadrangulations with two marked points at fixed distance. We then give 
in Section 3.2 the generating function for the average number of geodesies of length i 
in pointed quadrangulations with area n oo. In particular, for large i, we obtain an 
average number 3 x 2* of geodesies linking two given vertices at distance i. The average 
number of confiuent geodesies is computed in Section 3.3 and some of their refined prop- 
erties are discussed in Section 3.4, such as the average number of contacts between two 
weakly avoiding geodesies and the average area which they enclose. We finally discuss 
the case of strongly avoiding geodesies and show that, at large n, only a set of size 
oc n^/^ of exceptional points can be reached from a given vertex by two such strongly 
avoiding geodesies. In Section 4, we study a different scaling regime in which the size n 
of the quadrangulation and the length i of the geodesies tend simultaneously to infinity, 
keeping the ratio r = ijn}!^ fixed. This scaling is customary in the context of large 
quadrangulations and we recall in Section 4.1 the formula for the universal density p(r) 
of vertices at rescaled distance r from a given origin. We show in Section 4.2 that upon 
appropriate renormalizations, the average number of geodesies of rescaled length r, or 
of fc-tuples of confiuent geodesies is characterized by the same scaling function p(r). On 
the other hand, we show that the density of exceptional points is described by a new 
scaling function that we compute in Section 4.3. We finally discuss in Section 4.4 the 
area enclosed by two geodesies and show that this area is proportional to the total area 
n only in the case of strongly avoiding geodesies. We discuss our results and conclude 
in Section 5. 

2. Geodesies in quadrangulations and labeled trees 

2.1. Schaeffer's construction for planar quadrangulations 

In order to study maps with marked geodesies, we shall recourse to the well-known 
correspondence between rooted planar quadrangulations and so-called well-labeled trees 
||15|| in the form found by Schaeffer et al. The correspondence works as follows: 

starting from a rooted planar quadrangulation, i.e. a quadrangulation of the sphere 
with a distinguished oriented edge, we label each vertex by its graph distance from the 
origin, i.e. the vertex from which the root edge originates. Each vertex therefore receives 
a strictly positive integer label, except for the origin which receives the label 0. It is 
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(a) 



(b) 



Fig. 1: The two possible types of faces in a planar quadrangulation ac- 
cording to the distances of their incident vertices from a given origin: (a) 
simple faces and (b) confluent faces. To each face, we associate an edge of 
well-labeled tree (thick solid or dashed line). 

easily seen that the faces of the quadrangulation fall into two classes (see Fig.l): simple 
faces where the cyclic sequence of labels around the face is of the form z, i + 1, i + 2, 
i + 1 for some z > and confluent faces where it is of the form i, i + 1, z, i + 1. 




Fig. 2: An example (a) of rooted quadrangulation (thick or thin solid 
lines) and the associated edges (thick solid or dashed lines) in Schaeffer's 
construction. These edges form a planted well-labeled tree (b) . 

For each simple face, we select the unique edge of type i + 2 ^ i + 1 encountered 
clockwise around the face (see Fig.l-(a)). For each confluent face, we draw a new edge 
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linking the two vertices labeled i + 1 (see Fig.l-(b)). The procedure thus associates an 
edge to each face of the original quadrangulation and it can be shown that the graph 
made of these edges connects all the original vertices except the origin that remains 
isolated. Furthermore, it is planar, contains no loop and has one more vertex than 
edge, hence it is a plane tree (see Fig. 2 for an example). By construction, the labels 
on the vertices are strictly positive integers and vary by at most one along any edge 
of the tree. This property characterizes so-called well-labeled trees. The root edge of 
the quadrangulation selects a corner of the tree at its endpoint (labeled 1), at which we 
plant the tree. It was shown that the construction above provides a bijection between 
rooted planar quadrangulations and well-labeled trees planted at a corner labeled 1. 




(a) (b) 

Fig. 3: Construction of the quadrangulation associated with a given well- 
labeled tree. Corners are connected to their successors (a) by a set of 
non-crossing arches (dashed lines). These arches form the edges of the 
quadrangulation (b). 

The inverse construction can be realized as follows: starting from a well-labeled 
tree, we connect by an arch each corner with label i > 1 to its successor, that is the first 
corner labeled i — 1 encountered clockwise around the tree (see Fig. 3 for an example). 
Corners labeled 1 are then connected by arches to an extra vertex outside the tree. 
The fact that labels vary by at most one along an edge ensures that the arches can 
be drawn with no intersection, resulting in a planar quadrangulation while the labels 
simply record the graph distance from the added vertex. The arch from the corner 
labeled 1 at which the tree was planted to the extra vertex, when oriented backwards, 
defines the root edge of the quadrangulation. 
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2.2. Extension to quadrangulations with a marked geodesic 




(a) (b) (c) (d) (e) 

Fig. 4: The extension of Sciiaeffer's construction to quadrangulations with 
a marked geodesic. In (a), we only display the marked geodesic, while the 
underlying quadr angulation is symbolized by the grey background. We 
first unzip the geodesic path (b) and fill the created hole by squares (c). 
Applying the rules of Fig.l for the added squares selects all the edges of 
the left copy of the geodesic but one (d), creating a spine (thick line). 
The original faces of the quadr angulation give rise to well-labeled subtrees 
attached only to the left side of this spine, resulting in a spine tree (e) . 

We now turn to the study of quadrangulations with a marked geodesic, that is 
an oriented path connecting two vertices and having minimal length among all paths 
connecting these vertices. The vertex from which the geodesic originates will be referred 
to as the origin vertex. Clearly, when moving along a geodesic, the distance from the 
origin of the successive vertices increases by one at each step. Note that picking a 
geodesic of length 1 simply amounts to marking an oriented edge of the quadrangulation, 
in which case the object at hand is nothing but a rooted quadrangulation on which we 
may apply Schaeffer's construction. If the length i of the geodesic is larger or equal 
to 2, it is then convenient to adapt this construction to ensure that all the edges of 
the geodesic path but the first one are selected and form a path of length i — 1 on the 
tree. To this end, we decide to first "unzip" the geodesic so as to form a hole in the 
quadrangulation, bounded by two copies of the original geodesic (see Fig. 4), which we 
will refer to as the left and right copy upon viewing the hole with the origin at the 
top. We then fill up the hole with i — 1 new squares obtained by adding i — 2 parallel 
edges linking, for each k = — 2, the vertex of the right boundary at distance k 

from the origin to the vertex of the left boundary at distance k + 1. At this stage, it is 
important to note that both the unzipping process and the filling process do not affect 
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the distances from the origin on the quadrangulation. In particular, both copies of the 
original geodesic are actual geodesic paths in the modified quadrangulation and their 
vertex labels indeed correspond to the distance from the origin. We finally root the 
modified quadrangulation by selecting the edge — > 1 of the left copy of the geodesic. 

With this rooted quadrangulation at hand, we may now apply Schaeffer's construc- 
tion and transform it into a well-labeled tree. By construction, all the added squares are 
simple faces of which we select the lower most edge on the left boundary (see Fig.4-(d)). 
All the edges of the left boundary but the first one are therefore selected, thus forming 
a path of length z — 1 on the tree linking i vertices with increasing labels 1,2, ■ ■ 
We shall call this path the spine of the tree. On the other hand, none of the added 
edges is selected and the rest of the tree is therefore made of i (possibly empty) subtrees 
attached only to the left side of the spine, when viewing it as hanging from its first 
vertex. 

In view of this result, we are naturally led to define a spine tree as a particular type 
of well-labeled tree, made of a spine of length i — 1 joining i vertices with increasing 
labels 1, 2, • • • , i and of i well-labeled subtrees attached to the left side of this spine. The 
tree is planted at the first vertex of the spine and can be viewed as hanging from that 
vertex. 

With this definition, we see that applying Schaeffer's construction to our modified 
quadrangulation results in a spine tree. As we shall see below, not all spine trees however 
are obtained that way. This can be understood by first applying the inverse construction 
recalled in Sect. 2.1 to a spine tree obtained as above from a quadrangulation with a 
marked geodesic and following how the two copies of the geodesic path are reconstructed. 
For each vertex along the spine with label > 2, the associated corner lying just on the 
left of the incident spine edge pointing toward the root gives rise to an arch linking that 
corner to the preceding vertex on the spine (see Fig. 5). This set of arches, together with 
the arch linking the root corner to the extra origin vertex labeled 0, forms the left copy 
of the geodesic path on the modified quadrangulation. On the other hand, consider now 
the vertex labeled i at the end of the spine and its chain of successors (i.e. its successor, 
the successor of its successor, etc): this chain forms the right copy of the geodesic. 
Furthermore, one easily sees that, for j > 1, the corner labeled j + 1 lying on the right 
of the spine will be connected by an arch winding around the tree to the vertex labeled 
j of the chain. Completing the inverse construction, the resulting object is the desired 
quadrangulation where the two copies of the geodesic are now connected by z — 2 winding 
arches, with the rest of the quadrangulation in between. Upon continuous deformation 
on the sphere, we can recover the wanted zipper-like structure with the former winding 
arches now lying in between the two paths (see Fig. 5). This last structure is eventually 
zipped so as to recover the marked geodesic on the quadrangulation. 

Now we see that, in the last transformation, it is crucial that the two copies of the 
geodesic path obtained by the inverse construction have no common vertex (except for 
their first and last vertex). This is the case for a spine tree resulting from the modified 
Schaeffer's construction on a quadrangulation with a marked geodesic, but this is not 
true in general for an arbitrary spine tree. Indeed, for an arbitrary spine tree, it may 
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"^^S^^^-"^^ (b) 

Fig. 5: Starting from a spine tree obtained as in Fig. 4, we reconstruct 
the two copies of the geodesic path by linking corners to their successors 
(a). Corners on the left of the spine create the first copy (green short- 
dashed arches), while the chain of successors of the vertex i at the end of 
the spine creates the other copy (red long-dashed arches), the two copies 
being connected by winding (magenta dot-dashed) arches. They moreover 
surround all subtrees hanging from the spine hence, after completing the 
inverse construction, serve as boundaries for the rest of the quadrangula- 
tion (grey background in (b)). On the sphere, this configuration can be 
deformed into a quadrangulation with an internal zipper-like structure (c) 
which can be eventually squeezed into a marked geodesic (d). 

happen that the successor of a corner on the right of the spine lies itself on the spine 
(see Fig. 6), necessarily on its left side. If so, we will obtain two geodesic paths with a 
common vertex acting as a pinch point separating the rest of the quadrangulation into 
two parts (see Fig. 6). 

When applying the inverse construction to a general well-labeled spine tree with 
spine of length i — if we decide not to draw the winding arches, the most general object 
that we end up with is a quadrangulation with a geodesic boundary of length 2i, i.e. a 
boundary with a marked origin vertex and such that the two paths of length i joining 
this origin to its antipodal vertex are geodesic paths (see Fig. 7 for an illustration). This 
requirement is equivalent to demanding that the vertex antipodal to the origin be indeed 
at distance i from this origin in the whole quadrangulation, i.e. all paths in the bulk of 
the quadrangulation from the origin to that vertex have length larger or equal to i. The 
two geodesic paths forming the boundary may however meet at common vertices (with 
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(a) 



(b) 



Fig. 6: Appearance of pinch points in the inverse construction from a 
generic spine tree. If the successor of a corner labeled j + 1 on the right of 
the spine happens to be itself on the spine (a), this creates a pinch point 
at distance j from the origin in the associated quadr angulation (b). Here 
we have not represented the winding arches, as opposed to Fig. 5. 

necessarily the same label on both sides) or even along common edges. The construction 
above is clearly one-to-one. 

This constitutes bijection I between, on the one hand, quadrangulations with a 
geodesic boundary of length 2i (and possibly pinch points) and, on the other hand, 
planted well-labeled spine trees with spine of length i — 1. By convention, for z = 1, a 
well-labeled spine tree with spine of length is nothing but a well-labeled tree planted 
at a corner labeled 1. 

Moreover, there is a straightforward bijection 11 between quadrangulations with 
a geodesic boundary of length 2i having no pinch point and quadrangulations with a 
marked geodesic of length z, as obtained by gluing the two sides of the boundary around 
the sphere (see Fig. 7). 

We could easily characterize the subset of spine trees which are in correspondence 
with geodesic boundaries without pinch points. It is however much simpler to note that 
a quadrangulation whose boundary has k pinch points and length 2i may be uniquely 
decomposed into k -\- 1 irreducible components whose boundary lengths add up to 2i. 
Each of these irreducible components is either a single edge or is precisely a quadran- 
gulation with a geodesic boundary having no pinch point. This constitutes bijection III 
and is the final step that will allow us to relate the enumeration of quadrangulations 
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(a) 



(b) 



Fig. 7: A schematic picture (a) of a quadrangulation with a geodesic 
boundary of length 2i, as obtained from the inverse construction acting on 
a generic spine tree. The two geodesic paths bounding the quadrangulation 
may have pinch points in the form of common vertices and/or common 
edges. Quadrangulations with a geodesic boundary of length 2i having 
no pinch point (b) are in bijection with quadrangulations with a marked 
geodesic of length i. 

with a marked geodesic to that of spine trees. 

2. 3. Generating functions 

We now would like to count quadrangulations with marked geodesies. As usual, 
we will consider the generating function of such maps with a weight g per face. In 
the well-labeled tree language, this translates into a weight g per edge of the tree. We 
shall call Ri{g) the generating function for well-labeled trees planted at a corner labeled 
i. From Schaeffer's bijection, Ri{g) is nothing but the generating function of rooted 
planar quadrangulations or planar quadrangulations with a marked geodesic of length 
1. As mentioned in Ref.0, it can be shown that the quantity log (^Ri{g) / Ri-i{g)^ 
is the generating function for quadrangulations with two marked (and distinguished) 
vertices at distance i from each other (counted with symmetry factors). The generating 
functions Ri{g) satisfy the recursion relations 



for i > 1 with the convention Rq = 0. This relation simply states that a well-labeled 
tree planted at a vertex labeled i is formed by a collection of an arbitrary number of 
subtrees hanging from that vertex which are themselves well-labeled trees with a root 



R^{9) 



1 



(2.1) 



1 - g{R^_^{g) + R^{g) + R^+^{g)) 
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vertex labeled i — 1, i or i + 1. These recursion relations were solved in Ref.|ll], with 
the result: 



R,ig) = R{g)- 



where R{g) 



{1 - x{gy){l - x{gy+^) 
'l-xigy+^)il-xigy+^) 

6g 



(2.2) 



and x{g) 



1 - 24(7 - y/l - 12g + V6^/72g^ + 6g + ^1 - 12g - 1 

2 {6g + v/r^n^ - 1) 



In this expression, the functions R(g) and x{g) are explicit solutions of the equations: 



R{g) = l + 3g {Rig))' 



x{g) + 



x{,g) 



+ 1 



(2.3) 



All functions above have a singularity dX g = gcrit = 1/12, which governs the asymptotic 
behavior of the general term in their series expansions in g. From Ri{g), it is easy to 
obtain the generating function Zi{g) for well-labeled spine trees with a spine of length 
i — 1 as these trees are simply made of a spine with i attached well-labeled subtrees 
whose root vertices have labels increasing from 1 to i. We immediately get: 



Z.{g) = l[Rj{g) = R{gy 



{1 - xjgyjl - x{gy+^) 

{1 - x{gryi - x{gy+^) 



(2.4) 



For convenience, we decided not to assign the weight g to the edges of the spine. Indeed, 
in our construction, these edges correspond to the added squares which were not present 
in the original quadrangulation, and should therefore not be counted as real squares. 
From bijection I of previous section, Zi{g) is nothing but the generating function of 
quadrangulations with a geodesic boundary of length 2i, as defined above. 



(Q) 
(P 







Fig. 8: The set of all quadrangulations with a geodesic boundary of length 
4 and with or 1 face (in grey). 
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Fig. 9: A schematic picture of the configurations of quadrangulations with 
a geodesic boundary of length 2i as counted by Zi (generic case) or Ui (case 
with no pinch point), and of their relation (2.5). 

For instance, we have the expansion Z2 (g) = 1 + 5g + 32g^ + 234:g^ + . . . where the first 
two terms correspond to the configurations displayed in Fig. 8. 

Finally, we will denote by Ui{g) the generating function of quadrangulations with 
a geodesic boundary of length 2i made of two geodesic paths that have no pinch point. 
As we have seen, from bijection II, Ui{g) is precisely the desired generating function for 
quadrangulations with a marked geodesic of length i. From bijection III, the function 
Ui{g) is nothing but the irreducible part of Zi{g) obtained inductively via the relation 

U,{g) ^ Z,{g) ^ R,{g) 

^2 5) 

Ui{g) = Z,{g) -Y,Uj{g)Z,-j{g) , i>2 . 
i=i 

The subtraction term simply amounts to removing those quadrangulations for which the 
two geodesic paths have their closest pinch point to the origin at distance j, dividing the 
quadrangulation into a first irreducible part {Uj{g)) and a remaining part with possibly 
other pinch points {Zi-j{g)), as illustrated in Fig. 9. 

For instance, we have U2{g) = g + 10(7^ + 90g'^ + SlOg"^ + . . ., where the first two 
terms correspond to the configurations of Fig. 10 of quadrangulations with a geodesic 
boundary of length 4 having no pinch point, or equivalently, to the configurations of 
Fig. 11 of quadrangulations with a marked geodesic of length 2. 

2.4- Confluent geodesies 

With the above generating functions at hand, we have also access to generating 
functions for quadrangulations with several marked confiuent geodesies, namely quad- 
rangulations with say, k marked geodesic paths joining a common origin vertex to a 
common endpoint at distance i. We shall furthermore impose that these k geodesies 
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Fig. 10: The set of all quadr angulations with a geodesic boundary of 
length 4 having no pinch point, with one or two faces (in grey). 




Fig. 11: The set of all quadrangulations of the sphere with a marked 
geodesic (thick line) of length 2 and with one or two faces. These con- 
figurations are in one-to-one correspondence with those of Fig. 10. When 
drawing the quadrangulation in the plane, we choose as external face the 
one lying on the left of the — > 1 edge of the geodesic path. 

be weakly avoiding, in the sense that they do not cross each other, but still can have 
common vertices or even common edges. We shall also consider the case of strongly 
avoiding confiuent geodesies which are required to have no common vertex except for 
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their first and last endpoints. In both cases, we furthermore decide to distinguish one of 
the geodesies as the first one, which induces a natural linear ordering of the geodesies. 




(a) (b) (c) 

Fig. 12: Schematic picture of a configuration (a) with k weakly avoiding 
confluent geodesies (here k = 4). Upon unzipping (b) of the flrst geodesic 
path (here marked by an arrow) and continuous deformation on the sphere, 
we end up with a quadrangulation with a geodesic boundary (c) and k — 1 
geodesies in the bulk. It can be viewed as the juxtaposition of k quadran- 
gulations with boundary of the type depicted in Fig. 7- (a) with the global 
constraint that the two outermost boundaries cannot have common ver- 
tices. 



Now a quadrangulation with k such linearly ordered weakly avoiding geodesies may 
be transformed upon unzipping of the first geodesic and continuous deformation on the 
sphere into a quadrangulations with a geodesic boundary and k — 1 marked geodesies 
in the bulk of the quadrangulation (see Fig. 12). Note that these geodesies may very 
well have common vertices or edges among themselves and with the boundary if the 
corresponding geodesies had such common vertices or edges on the quadrangulation. 
Nevertheless, there is a global requirement that the two boundaries themselves cannot 
have common vertices. Without this requirement, the generating function would simply 
be [Zi{g)j by viewing the configuration as simply made of k glued objects counted by 
Zi{g). To account for the global constraint, we again simply have to extract the irre- 
ducible part of (Zi((7))'^, leading to a generating function ul^\g) for quadrangulations 
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with k marked weakly avoiding confluent geodesies of length i given by 



Ui'\g) = {Z,{g))' 

j=i 



Finally, the situation is even simpler in the case of k strongly avoiding geodesies 
since in this case, upon unzipping all geodesies, we obtain exactly k quadrangulations 
with a geodesic boundary with no pinch point. The generating function uj;''\g) for 
quadrangulations with k marked strongly avoiding geodesies of length i is therefore 
given by 

U^''\9) = {U.{9)t ■ (2-7) 



3. Statistics of geodesies in large quadrangulations 

3.1. From combinatorics to statistics 

Let us now see how to use the above generating functions to study the statistics of 
geodesies in quadrangulations. We shall be interested in averaging over all quadrangu- 
lations with a fixed number n of faces. As customary, our enumeration results translate 
naturally into averages for an ensemble where every quadrangulation is drawn with a 
probability proportional to its inverse symmetry factor, i.e. to the inverse of the order of 
its automorphism group. In the limit n ^ oo of large quadrangulations, it is expected 
that the symmetry factors become irrelevant, so that the asymptotic statistics is the 
same as for the uniform probability distribution. 




1 1/2 1 1/2 1/2 1 1 1/2 

Fig. 13: The set of all pointed quadrangulations of the sphere with one 
and two faces, where we have indicated the marked vertex by a label 0. 
For each configuration, we indicated its inverse symmetry factor. These 
add up to = 3/2 for 1 face and to ^2 = 9/2 for 2 faces. 



In our case, a natural quantity to look at is the number of geodesies emerging 
from a given origin. In this context, it is natural to consider the ensemble of pointed 
quadrangulations with n faces, i.e. quadrangulations with a marked origin vertex, again 
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weighted by their inverse symmetry factor. The probabihty of a given pointed quad- 
rangulation is then equal to 1/(32^) where S is the symmetry factor and is the 
partition function, which reads 

3. (2n) 



2. 



2nn + 1 



(3.1) 



For instance for n = 1 and 2, the results Zi = 3/2 and Z2 = 9/2 are illustrated in 
Fig. 13. Note that, in contrast with rooted quadrangulations, pointed quadrangulations 
may have non-trivial symmetries, hence is not necessarily an integer. 



o 

T 
O 

1x1 1/2x0 





1/2x0 





1x2 













(2) J 





1x2 




Fig. 14: The average number of geodesies of length 2 emerging from the 
origin of pointed quadrangulations with 1 and 2 faces. This average is 
obtained as the sum over pointed quadrangulations of the actual number 
of geodesies (indicated in green on the right) times the symmetry factor of 
Fig. 13 (indicated in red on the left), leading to 1 (for n = 1) and 10 (for 
n = 2), and dividing by the partition function. 
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In this ensemble, we can interpret the quantity Ui{g)\gri/ as the average number 
of geodesies of length i emerging from the origin. For i = 2, we have U2{g) = g + ^Og^ + 
. . ., so that, for n = 1 and 2, this average number reads respectively 1/(3/2) = 2/3 and 
10/(9/2) = 20/9 (see Fig. 14 for an illustration). 

Alternatively, we can be interested in the ensemble of quadrangulations with 
two marked vertices at distance i from each other, again weighted by symmetry 
factors. In that mentioned above, the partition function is known to be 

log (^Ri{g) / Ri-i{g)^\gn . In this ensemble, the average number of geodesies joining the 
two marked points is simply Ui{g)\gri divided by this different normalization. 

3.2. Average number of geodesies in large quadrangulations 

We are now interested in the large n limit of the statistics of geodesies. To this 
end, we first need to extract the coefficient of g"" in the various generating functions of 
previous section. For an arbitrary series X(g), this coefficient is obtained via a contour 
integral around 0: 

X{g)\g. = ±-j^^X{g). (3.2) 

In the case of our generating functions, the large n behavior is governed by their singu- 
larity at the critical value gcrit = 1/12. To capture this singularity, it is convenient to 
first change variable from gloV = gR{g), with R{g) as in Eq.( p.2|) , such that 

g = g{V) = V-3V^ , (3.3) 

leading to 

1 r dV 



X(g)\„^ = (p -j(l-6V)X(g(V)) . (3.4) 

^^^'^ 2i7r7 (V(l-3y))"+' ^ ^ 

At large n, the dominant behavior is given by the vicinity of saddle point at Vcrit = 1/6. 
Writing 

V = l(^ + ^-D : (3-5) 



6 V nV2^ 

the contour integral over V becomes at large n an integral over real values of ^: 

(..^.0(l))x(l(.4)) . (3.0, 
For X = Zi or Ui, we will find an expansion of the form 

with real coefficients A, C and and with a vanishing ^jn}!'^ term. We have introduced 
the prefactor 2/3 for later convenience. The even powers of ^ in (|3.7|) do not contribute 
to (|3.6|) by parity and we get the leading behavior 

12'^ /"+°° 2 2 12" 
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As explained in previous section, whenever X{g) enumerates objects emerging from an 
origin, the corresponding average number (x) in the ensemble of pointed quadrangula- 
tions is obtained by dividing the above result by: 

12" 

so that the asymptotic average simply reads 

{x) = D . (3.10) 
Using the explicit form (|2.4|) for Z^, we find immediately the expansion 



Z, = A, + a^ + -iD,Arj^ + .. 

-.1. . 2^(^ + 3) 
with Ai - — 



D, 



3(z + l) 
2H{i + 3){f + 6i + 3) 
30(2 + 1) 

2H{i + 2){i + 3){i + 4:){3f + 12i + 13) 



(3.11) 



420(z + 1) 

from which we obtain the large n asymptotic behavior 

7(M 12- 2-z(z + 2)(z + 3)(i + 4)(3z^ + 12z + 13) 

- 7^ 840(z + 1) ^^-^^^ 

As for Ui, we have a similar expansion 

U, = a,+-f,^ + -i5,^ + ... (3.13) 
n 3 n"^/^ 



where the coefficient Si gives the large n behavior of Ui via the general relation (|3.8| ). 
In particular, from ( |3.10| ), the average number of geodesic paths of length i emerging 
from the origin vertex in the ensemble of pointed quadrangulations of fixed large size is 
directly given by Si. The coefficients ct^, 7^ and Si can be related to Ai, Ci and Di via 
the relation ( |2.5| ), namely 

i-l 



Ai -"^ajAi 

i-l 

l^ = a-J2 + IjA^-j) (3.14) 

i-l 
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These relations are better expressed by introducing generating functions for these coef- 
ficients: 



A{t) = J2 C{t) = J2 C^t\ D{t) = J2 Drf 

i=l i=l 1=1 

oo oo oo 

(f)^5^a,f, ^{t)^Y.^,t\ 5(t)^5^5,f 



(3.15) 



a 

i=l i=l i=l 

as they translate into: 

Mt) , . Cit) . D(t) 

l + A{t) (1+A(t))2 (1+A(t))2 

From the explicit values of Ai, Ci and Di, we can get explicit values for all these 
generating functions, leading to: 



m = 
m = 



t{6t - 2) - (1 - 2t) log(l - 2t) 

t - (1 - 2t) log(l - 2t) 
3t {2t {20t^ - 34t2 + m - l) - (1 - 2t)4log(l - 2t)) 

5(1 - 2t)2 (t - (1 - 2t) log(l - 2t)f 
3t (2t (3 + 177t - 412^2 + 708t^ - 624^^ + 224t^) + 3 (1 - 2tf log(l - 2t)^ 

70 (1 - 2t)^ (t - (1 - 2t) log(l - 2t)f 

80 9 o 366208 , 3998176 , , 

4t + —t^ + 132^3 + + + O (t^) . 

3 675 2025 ^ ^ 



(3.17) 

The function S{t) is the generating function for the average number Si of geodesies of 
length i emerging from the origin vertex in a large pointed quadrangulation. From the 
last line of (|3.17|) , we directly read of the values of this average number for a length 



i = l,2,---,5. 

We can easily estimate the large i behavior of ai, 7^ and 6i by looking at the 
singularity of a{t), ^{t) and d{t) for t 1/2. Writing t = 1/2(1 — 77), we have the small 
Tj behavior 

a Q(l -V)^ = 1 - 377 - 6r7^1og(?7) + . . . 



2' "J 7774 

Note that this last behavior follows immediately from the small 77 behaviors A{t) 
1/(377) and D{t) ~ 6/{7r]^) and from the relation (|3.16|) . From ( |3.18| ), we deduce: 



i — ^00 12 i — ^00 (-.7 9 , , i — ^00 ^i 9,0 /f-, ^ ^\ 

a, ~ 2* — , 7i ~ 2'- I, di ~ 2* - . (3.19) 

1 O I 
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Prom the discussion of previous section, to go to the ensemble of quadrangulations with 
two marked vertices at distance z, we simply have to gauge Si by the average number of 
points at distance i from the origin in large pointed quadrangulations. This number is 
given by {3/7)i^ for large i (see for instance Eq.(5.2) in Ref.p), hence we deduce that 
the average number {Q)i of geodesies joining two fixed points at distance i from each 
other is simply given by 

(G), -3x2* (3.20) 

at large large i. The leading exponential i-dependence can be heuristically understood 
as follows: note that choosing a particular geodesic leading to a given vertex simply 
amounts, when going backwards from that vertex, to choose recursively, for each vertex 
at distance i a predecessor vertex at distance z — 1, until we eventually reach distance 
0. Now one can easily show that, for a vertex at distance i ^ 1 from the origin of a 
random quadrangulation, its average number of neighbors at distance i — 1 is 2. There 
are thus on average 2 choices for each of the i steps leading back to the origin, which 
explains the 2* factor by neglecting correlations. The prefactor 3 precisely accounts for 
these correlations as well as for the finite i corrections to the asymptotic number of 
choices 2. 



3. 3. Average number of confluent geodesies 

We can as well obtain the average number of A;-tuples of (linearly ordered) weakly 
avoiding confluent geodesies of length i by considering the expansion: 

(zo'^ = Af) + cf)^ + ^iDf 4^ + ... 

n 3 n'^/^ 

(3.21) 



with Af^ 


= {A.t 






= k {Ai] 






= k (A,; 





where the coefficients Ai, Ci and Di are given by Eq.( |3.11| ). In particular, the generating 
functions A^^\t) = ^-^^ , C^''\t) = E,>i cf and D^'^^t) = J2i>i ofh' now 

have a singularity at t = 1/2^^. 



(k) 

Similarly, the function t/- has the expansion 



+l^ - + ^ + ^3.22) 

where, as before, the coefficients oii'\ '^[^^ and 5^^^ are determined through their gen- 
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erating functions: 

i>l 



2 



(3.23) 



2 



Writing t = 1/2^^ (1 — 77) and expanding at small ry, we have now the dominant be- 
haviors A^^\t) ~ 1/(3^=77) and D^^\t) ~ 6/c/(3'=-^ 7?7^), leading directly to S^^\t) ~ 
(fc3^~^ 54)/(777^). From this singularity, we immediately deduce the large i behavior of 

'a~(2')'*:3'-'^i3 (3.24) 

which gives the average number of fc-tuples of (linearly ordered) weakly-avoiding con- 
fluent geodesies of length i emerging from the origin of a large pointed quadrangulation. 
Again, when compared with the average number (3/7)i^ of vertices at distance i from 
this origin, this gives an average number 

{gk)^ ~ A;(3 X 2^ (3.25) 

of such /c-tuples between the two marked points in the ensemble of quadrangulations 
with a large fixed size and two marked points at distance i from each other. 

Let us make a few remarks on this particularly simple result. First, the leading 
exponential z-dependence 2*^^ is surprisingly the same as for k independent geodesies. 
One could have expected a reduction of the entropy factor due to the constraint of non- 
intersection between geodesies. Indeed, as we shall see later, weakly avoiding confluent 
geodesies have a number of contacts which is proportional to their length. This apparent 
paradox can be explained as follows: as before, we can estimate the number of geodesies 
by looking at them backwards and choosing recursively for each vertex at distance i one 
of its predecessors at distance z — 1. For a vertex visited by a single geodesic, this 
gives 2 choices on average as before. When a number £ of geodesies come in contact at 
a vertex with p predecessors, the number of choices is now (instead of for 

freely intersecting geodesies). Now the probability of having p predecessors for a vertex 
at a large distance from the origin can be easily shown to be 1/2^ (with p varying from 
1 to 00), leading to an average value ((^^^~^)) = 2^, i.e. the same contribution as if 
the i geodesies were visiting independent vertices. In other words, the reduction of 
choices among predecessors due to non-intersection is exactly compensated by the effect 
of correlations for the number of predecessors at a given vertex. 

Even more surprising, the sub-leading prefactor 3 in Eq.( p.2Ci|) is simply raised to 
its k-th power as if the effect of small i corrections and correlations along the geodesies 
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could be treated independently for each geodesic path. Finally, we have in Eq.( p.25|) a 
global prefactor k which can be understood as a spontaneous symmetry breaking effect 
as follows: the k confluent geodesies are linearly ordered and any pair of neighboring 
geodesies encloses a particular domain of the quadrangulation, with k such domains. 
As we shall see, for quadrangulations with a large number n of faces, only one of these 
k domains contains most of the area of the quadrangulation, i.e. has a number of faces 
close to 71, while all the other domains have an area negligible with respect to n at large 
n. There are precisely k choices for the domain in which most the area will lie, hence 
the global prefactor k. 

3.4- Other statistical properties of confluent geodesies 

By slight refinements of the above arguments, we have access to more involved 
statistical properties of confluent geodesies. For instance, we may compute the asymp- 
totic (large n) average number {c)i of contacts between two weakly avoiding confluent 
geodesies of length z, with result: 

{ch ~ ^ (3.26) 

at large i. A proof of this result is presented in Appendix A using the continuous 
formalism of Section 4. As announced, this number is proportional to the length of the 
geodesies which indicates that the two geodesies "stick" to one another. 




(a) (b) (c) 

Fig. 15: Schematic picture of a configuration (a) with two weakly avoiding 
confluent geodesies. The two geodesies delimit two (possibly disconnected) 
complementary domains (here in dark and light) lying between them. The 
same configuration (c) obtained upon unzipping of one geodesic path (b) 
and continuous deformation on the sphere. 

Another quantity of interest is the area (number of faces) enclosed by these two 
geodesies. More precisely, let ni and n2 be the areas of the two domains delimited by 
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the geodesies in the quadrangulation, with ni + n2 = n (see Fig. 15 for an iUustration) . 
By symmetry, we have (ni) = (n2) = n/2 and a good measure of the dispatching of the 
area is given by the correlation (nin2). For two geodesies of length i, we find: 

{nin2)i ~ nOi (3.27) 

for n — s> oo and i finite, with a finite coefficient 9i. This implies a very dissymmetric 
dispatching where the smaller domain has an average area of order 9i only while almost 
all of the area lies in the larger domain. For large i, we obtained the leading behavior: 

27 

0. ~ ^/ . (3.28) 

A derivation of this result is given in Appendix B by use of the continuous framework of 
Section 4. The above dissymmetry generalizes easily to the case of domains delimited 
by k confluent geodesies. We find that only one of these domains has an extensive area, 
meaning that the geodesies remain close to each other. 

In view of the above results, we expect very different statistics for strongly avoiding 
confluent geodesies where contacts are forbidden. The asymptotic number of /c-tuples 

~ (k) 

of such geodesies of length i is obtained from the expansion of U^^ , namely: 

with ^ = {aif ^^^g) 
= k (a.)'"' ■ 



From the large i behavior ( |3.19| ), we get 



consisting of the same degeneracy factor k{3 x 2*)^^ as in Eq. (|3.25|) multiplied by a factor 



(3 4^~^/7) which we interpret as counting the average number of vertices that are 

indeed reachable by k strongly avoiding geodesies. For /c > 2, this factor goes to zero at 
large i, which indicates that distant vertices cannot be reached by more than 2 strongly 
avoiding geodesies. For = 2, we have for each (large) value i of the distance, an average 
number 12/7 of vertices reachable from the origin by 2 strongly avoiding geodesies. We 
shall call such vertices exceptional points. 

Let us now consider the areas ni and 77-2 of the two domains delimited by 2 strongly 
avoiding confluent geodesies of length i. Again we want to compute the correlation 



(nin2). This is done by considering the generating function {9-j^Ui{g)j as the action 



2 
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oi g-j^ amounts precisely to giving a weight proportional to the area. When going to 
the ^ variable, we have at leading order in n that 9^ = hence, from the expansion 
( p.l3|) , we get the expansion 

9YgU^i9)) =7? + 2i7.'J.^ + ... (3.31) 



Performing the integral over ^, the first term vanishes by parity and we get the dominant 
behavior 

d \^ 12" 

This is to be compared with the asymptotic behavior 

iU^{g)f ~ TT^^^f ^ = -^^^^^^ ■ (3-33) 

Taking the ratio of (|3.32| ) and (|3.33| ), we get the asymptotic average 

{nin2)i = n— ~ —m . (3.34) 
ai 2(j 

Again, almost all of the area is carried by a single domain, while the smaller domain now 
has an area of order i^. This result is valid in the limit n ^ 00 and i finite but large. 
It will be extended in the next section to the scaling limit of large quadrangulations, 
which is obtained by taking z oc n^/^. In this limit, we have oc n, which means that 
the smaller area is of the same order as the larger one, i.e. the two domains contain a 
finite fraction of the total area. 



4. Continuum limit 

4.1. Reminders 

So far the results which we presented were obtained in the limit of a size n of 
the quadrangulations tending to infinity, keeping a finite value for the length i of the 
geodesies. On the other hand, it is well known that an interesting scaling limit may be 
obtained by letting both n and i tend to infinity, while keeping the ratio i/n^^'^ fixed. 

As a reminder, let us first consider the average number of vertices within a ball 
of radius i = rn^^^ centered at a given origin in a random quadrangulation of size n. 
In the limit n 00, this average number behaves as n$(r), with $(r) a finite scaling 



function. This function was computed in Ref.[|TT|, with the result: 







(cosh(ri/3|') — cos{ry/3^)) 
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(a) 



(b) 



Fig. 16: The distribution function $(r) measuring the probability that 
a vertex be at a rescaled distance less than r in the scaling limit of large 
quadrangulation. Its derivative p(r) is the density of vertices at rescaled 
distance r. 

It increases from $(0) =0 to $(+oo) = 1 as it should and its derivative 

p(r) = $'(r) 



(4.2) 



may be interpreted as the average density of vertices at a rescaled distance r from a 
given origin. Both functions $(r) and p(r) are plotted in Fig. 16 for illustration. We 
have in particular the asymptotic behaviors: 



p{r] 



3 
7 



e-3(3/8)^ 



(4.3) 



The rest of this section is devoted to the computation of similar scaling functions char- 
acterizing the statistics of geodesies in the same scaling limit. 



4-2. Distribution of geodesies of length r 

Let us first consider the case of quadrangulations of size n, with a single marked 
geodesic of length z, as counted by Ui{g)\gn. To estimate this quantity in the scaling 
limit, we first evaluate Zi{g)\gn. from its explicit form (|2.4|) . For i = rn^^^, we have the 
expansion 



2* V 12 V n 



3 7iV4 '^^\n^/2 



with J^{r,^) 



- - i £ coth r 
3 ^ ' 




(4.4) 



25 



which, upon substituting in (|3.6|), gives the asymptotic behavior 



2^ 



n — »-oo 



g„ Trn*^/* Jo \ 3 V 3 cosh(rv34) - cos(rv34) / 



In particular, for r ^ 0, the integral behaves as leading to 

12"/(0r7i^/'^) X (r^/280), consistent with Eq.(|?T|) in the regime K z -C n^/^. 
As before, we may obtain IJi from via equation (|2.5|) , which we rewrite as 



lJ{t) = (4.6) 
1 + Z{t) 

upon introducing the generating functions 

oo oo 

= = (4.7) 

1=1 1=1 

with a weight t per geodesic step. The scaling limit i oc n^^^ oo is now captured by 
taking 

t = ^e-^""'^\ (4.8) 

In this regime, we cannot simply plug the asymptotic form (|4.4| ) of Zi in Z{t) as J^(r, ^) 
has a pole at r = 0. A simple way to work around this problem is to subtract from Zi 
its critical value ^^(1/12) = Ai, as defined in Section 3.2, and to write 

Noting that the right hand side has no more pole, the sum over i can be approximated 
by an integral over r, leading to: 

Z(t)-i(t) = y^°°dre-'-(^^(r,0-^)+o(^) , (4.10) 

with, as before, A{t) = From the explicit form ( |3.11| ), we find the explicit 

expansion of A{t): 

Mt) = ^--^og(^)-'- + o('] . (4.11) 
Plugging these expansions into ( [4.6| ) , we deduce the expansion of tj if) 



+ . . . 

(4.12) 
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where we have identified the expansion of a{t), as defined in Section 3.2, and with 
JF" = d'^T jdr'^ . Going back to Ui{g), this gives an expansion 

^(f^^-«^) = ^(9-^>'^)-^)+^Q) ' (4-13) 
which, from the asymptotics ~ of Eq. ( |3.19| ), yields the scahng form 



n 



Surprisingly, when comparing with Eq.( |4.4| ), we see that, in the scaling limit, extracting 
the irreducible part simply amounts to differentiating twice with respect to r. Upon 
substitution of ( [4.14 ) in Eq.( |3.6D , we obtain the desired asymptotic scaling behavior 



El 

2* 



^^7/4 Qj.2 1 3 Y 3 cosh(rv/3^) — cos{r^/^) 

12^ 12" 

X 3$'(r) = X 3p(r) 



20Fn7/4 ' ' 2v/^nV4 

(4.15) 

Up to a factor of 3, this last expression is nothing but the number of quadrangulations 
with simply a marked vertex at distance i = rn}l^ from a given origin. We therefore 
recover in the scaling limit the same ratio 3 x 2* as obtained in Eq. (|3.20|) for 1 -C z -C n^/^ 
between geodesies of length i and vertices at distance i. Up to this degeneracy factor, 
the density of geodesies of rescaled length r coincides with the density of vertices at 
distance r. 

We can easily play the same game with /c-tuples of weakly confluent geodesies. 
Using the expansion 

^(Zf - - {nr.O -l)^0 (-i,) (4,16) 

— 1/4 , 

we now get, for t = e~^'^ /2 , 



leading to 

12^ 



(2 



n — >oo 



2V^n7/4 



A;3VW ■ (4-18) 



Again, in the scaling limit, the density of /c-tuples of weakly avoiding confluent geodesies 
of length r coincides with the density p{r) upon a simple renormalization by a degeneracy 
factor /c(3 x 2*)^^, in agreement with Eq.( |3.25| ). 
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In the case k = 2 and by a slight refinement of our generating functions, we may 
compute the average number of contacts between two weakly avoiding confluent geodesic 
of rescaled length r. This calculation is detailed in Appendix A, with the result: 

{c)r = n^l^ I , (4.19) 



which extends the flnite i linear dependence displayed in Eq. (|3.26| ) to the whole range 
of r in the scaling regime. 

4-3. Exceptional points 

The situation is more interesting for exceptional points as it will lead us to a new 
scaling function. We consider again the case of two strongly avoiding confluent geodesies 
of length i, as counted by Uf. From Eg. (14.14]) , we now get for i = rn^^^ the scaling 
behavior 



1 on f+oo 

/ d^^e-«^9^"(r,0)' 

^ — OO 



(4.20) 



(2 X 3^) a(r) 



where 

a(r) = ^/ rfee'e-^-<'8^f-— ^ + 2 



432 r [g (^g-l)^ I , (c5-l)(cc~-4) 1 

^io I (c-cT (c-cT (c-5)2J (4.21) 

with c = cosh(rA/3^) , c = cos(rA/3^) . 




Fig. 17: The scaling function a(r) counting the number of exceptional 
points at a rescaled distance r form the origin in the scaling limit of large 
quadr angulations. There are on average n}l'^a{r)dr such exceptional points 
at a rescaled distance between r and r ^ dr. 
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We have factored out in Eq.( 4.20| ) the asymptotic number 12"'/(2-y7rn^/^) of quad- 
rangulations with a given origin as weU as a degeneracy factor 2x3^ for two confluent 
geodesies. We may then interpret the remaining finite scahng function cr(r) as counting 
the number of exceptional points at a rescaled distance r from a given origin. This 
function is plotted in Fig. 17 and has the following limiting behaviors: 

a(r) 12/7 

(4.22) 

a{r) g-3(3/2)2/3^4/3 



Note that the small r limit corroborates the finite i result of Section 3.4. Finally, 
upon summing over i, the total number of exceptional points is simply given by 
n"*^/* /o°° dra{r). 

4.4- Area between two geodesies 

Let us finally address again the question of the dispatching of the area between 
geodesies, now in the scaling limit. Considering two non-intersecting geodesies, we shall 
call LOi and uj2 the fraction of the total area n carried by the two domains which they 
separate, with of course UJ1+UJ2 = I- A convenient way to access the average properties 
of iVi and lv2 is to introduce extra weight parameters Ai and A2 coupled to the area of 
the two domains. In the case of weakly avoiding geodesies, we therefore consider the 
refined generating function Zi{\ig)Zi{\2g) counting quadrangulations with a weight 
Ai^r (resp. \2g) per square in the first (resp. second) domain. Now we may explore the 
scaling limit of Zi{Xig)Zi{X2g) by taking Ai = 1 + (9(l/n), namely by writing 

^ = 77^ H 



12 V n 



with iii = ^i(Ai,^) such that /Ui(l,^) = ^. At large n and for i — rv}!'^^ we have the 
expansion: 

-^Z,{Xig)Z,{X2g) = + — T7T + ^ rTT^ ■ (4-24) 



(202-'' '^'^' "^' 32' 3^1/4 ' ^ \n^/^ 

he corr 
with expansion 



Taking the corresponding irreducible part leads to a generating function {g; Ai, A2) 



1 rr(2). . 27(r-(r,^i)+rXr,^2)) , f I \ 
j-^U^ {g;X^A2) = + 0^-] (4.25) 
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and we end up with 



(2 



i\2 



n — >oo 



12" 

774 



irn 



27 / dCCe-^ {r'{r,^^^)+r'{r,^^2)} ■ (4.26) 
J —00 



To extract {ui) {i = 1,2), we must apply the operator {Xi/n){d/dXi) and let 1, 
i.e. Hi ^. At leading order in n, we have {Xi/n){d/dXi) = l/{2ni){d/d/j,i) and we 
thus get 



('^l)^ = Mr = "'Zoo,.. :L. \/ = o (4-27) 



as seen by an integration by part. This average result simply expresses the symmetry 
between the two domains and does not tell us anything about the actual dispatching of 
the area. As before, this dispatching is better measured by the correlation (t<;iCi;2)- This 
latter average requires the action of the operator l/{4fxiH2)d'^/dfXidfX2 on the integrand 
in ( [1.26|) , which produces a vanishing result, namely 

{uiU2)r = (4.28) 

in the limit n ^ 00. In other word, the fraction of the area of the smaller domain tends 
to at large n. A closer look at the corrections to scaling leads to {uJiuj2)r = C(l/n^/^), 
with an r-dependent coefficient. An integral formula for this coefficient is presented in 
Appendix B. In particular, for small r, it behaves as 27r'^/100, which is precisely the 
result announced in Eq.( p.28|) when i = rn}l^. 

This situation is to be contrasted with the case of two strongly avoiding geodesies, 
for which we must consider directly the expansion of Ui{Xig)Ui{X2g), namely 



1 



—U,{Xig)U,iX2g) 



IT 



/ di i e-^ 81^"(r, Hi)r'{r, ^2) • (4.29) 
J —00 



Applying again l/{2fXi)d/dfXi, we recover (^1)7. = (^2)7- = 1/2, while we now obtain 



2 



A(r) 



(uiUJ2)r = T ^ TT- = ^4 (4.30) 



where 



V^Jo ^ V dU V (c-c)6 (c-c)4 



{c-cy 

_ (cc-1)^ (cc- l)(2cc-5) cc-4 

' 'c-c)S (c-c)*' (c-c) 

(4.31) 
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(a) (b) 

Fig. 18: The scaling function A(r) of Eq.(4.31) and the ratio \{r)/(j{r) 
measuring the correlation {uJiLjj2)r for the fractions of area wi and ljJ2 of the 
two domains delimited by two strongly avoiding confluent geodesies. The 
average is for geodesies leading to a given exceptional point at distance r 
from the origin. For small r, \{r)/a{r) ~ (3/20) r"^, which means that the 
larger fraction of area is of order 1 and the smaller of order r^. For large 
r, \{r)/a{r) 1/4, which means that both fractions tend to 1/2. 



with c = cosh(rA/3^), c = cos{ry/3^), s = sinh(ry3^) and s = sm{ry/3^). This function 
is plotted in Fig. 18 together with the average value (ci;it<;2)r = X{r)/a{r). At small 



r, the quantity between curly brackets in (|4.31|) behaves as (36/35)^ r , leading to 
A(r) ~ (9/35) and {ujiu!2)r ~ (3/20) r^. We thus recover an area of order n in 
the larger domain and an area of order (3/20) nr^ = (3/20) in the smaller one, in 



agreement with the finite i result (|3.34|) . For an arbitrary finite value of r, both uJi and 



L1J2 are finite, leading to two domains whose area is proportional to n. This indicates 
that exceptional points are reachable by two geodesies which are truly distinct, even in 
the continuum limit, as they separate macroscopic domains of extensive area. Finally, 
for large r, A(r) vanishes like e"^*^'^/^-* ^^"^ as cr(r), with a limiting ratio 



A(r) ,.^00 1 
cr(r) 4 



(4.32) 



Now we note that we have necessarily uj\uj2 < 1/4 since uj\ and uj2 are real positive 
numbers adding to 1 and we have equality only for cui = c<;2 = 1/2. The large r limiting 
average value (^1^2) = 1/4 therefore indicates that the exceptional points sitting at a 
large rescaled distance r from the origin separate the quadrangulation in two domains 
of equal area n/2. Note that such exceptional points may appear only in pointed 
quadrangulations with a rescaled radius larger than r, which are exponentially rare. 
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5. Conclusion 

In this paper, we have presented a number of exact results on the statistics of 
geodesies in quadrangulations of large size n. At a discrete level, the number of geodesies 
between any two vertices at a distance i from each other grows like 3 x 2* at large i. 
For a generic pair of vertices, any two geodesies linking these vertices stick to each 
other in the sense that their number of contacts is proportional to i. Moreover, they 
separate the quadrangulation in two very asymmetric domains, with the larger domain 
containing most of the area of the quadrangulation while the smaller one has an area of 
order <^ n. This symmetry breaking extends to the case of fc-tuples of geodesies, for 
which only one of the k domains which they form in the quadrangulation has an area 
of order n. 

Upon factoring out the entropy factor 3x2*, a sensible scaling limit can be reached 
in the regime i = rn^^"^. In this limit, the density of geodesies of rescaled length r 
simply coincides with the density of vertices at distance r. Similarly, for /c-tuples of 
geodesies, we obtain the same density upon dividing by the entropy factor /c (3 x 2')^, 
where the first k prefactor is an expression of the above-mentioned symmetry breaking 
phenomenon. A possible interpretation of these results is that there is essentially a 
unique geodesic path between any two points in the scaling limit and that the entropy 
factors simply account for a degeneracy at a microscopic level. More studies are however 
necessary to validate this picture. In particular, we noted that the area enclosed by two 
geodesies, although negligible with respect to n, still scales as n"^/^, i.e. faster than the 
length i oc n^/'^ of the geodesies. Viewing the two geodesies as forming a sequence of 
blobs without contact, this means that at least one of these blobs has an area larger 
than n^/^. It is not clear whether the observed scaling is sufficient to guarantee that we 
can simply ignore the domain trapped between the geodesies in the scaling limit. 

In addition to this generic behavior, we also found that, for a given origin vertex, 
there is a number of order n-*^/^ of exceptional points that can be linked to this origin 
by two truly distinct geodesies with no contacts on the way. Both domains separated 
by the geodesies have an area proportional to n in this case. We explicited the density 
of these exceptional endpoints in the scaling limit as a function of the rescaled distance 
r but a more detailed analysis of their correlations would be desirable to understand 
better the geometry of this set of points. 

It would be nice to extend our results to other families of maps. We expect in 
general that the precise value of the degeneracy factors counting discrete geodesic paths 
will depend on the family at hand but that all our results in the scaling limit should 
remain unchanged. This holds for maps in the universality class of so-called "pure 
gravity". Studying models of maps with matter degrees of freedom is another story 
which requires new insights. 
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Appendix A. Average number of contacts between two confluent geodesies 

We compute here in the scahng hmit i = rrt}!^ the average number of contact 
vertices between two weakly avoiding confluent geodesies of length i. To this end, it 
is convenient to introduce the generating function Zf^{g) of quadrangulations with a 
geodesic boundary of length i and with a marked pinch point. We have the relation: 

Zf = J2z,Z,_,, (A.l) 

which simply expresses that the marked pinch point, at arbitrary distance j, separates 
the quadrangulation with a boundary in two parts of length j and i — j- 

To treat the case of two geodesic paths with a marked contact, we simply have to 
consider the generating function 2ZiZf^, which counts quadrangulations with a geodesic 
boundary of length 2z, with a marked geodesic path in-between, and with a marked 
contact between this geodesic path and either side of the boundary. Note that if the 
two boundaries come into contact, the corresponding pinch point is actually counted 
twice. For a proper enumeration of quadrangulations with marked geodesies, we must 
as usual take an irreducible part. In the same way as the irreducible part uj^"^^ of Zf 
counts quadrangulations with two weakly avoiding confluent geodesies, the irreducible 
part u'f'^^^ of 2ZiZ\^ counts these same objects with a marked contact. It is obtained by 
removing conflgurations where the two boundaries come into contact, which translates 
into the relation 

i-1 

Vf^^^ = {2Z,ZfP) - {U^^^ X 2Z,_,Zf^^. + t/j') X 2Zl^ + uf^" x Z^,) , (A.2) 

j=i 

where the subtracted terms correspond to conflgurations where the two boundaries have 
a flrst contact at a distance j from the origin. This contact may be before the marked 
pinch point (flrst term in the sum), at the marked pinch point (second term) or after 
the marked pinch point (third term). 

As before, the above relations transform into simple relations for the generating 
functions 



z(t) = zfi , z(2) (t) = Y zft" , zpp{t) = Y 

i=l i=l i=l 

oo oo oo 

i=l i=l i=l 

Indeed, Eq.( |A.l| ) translates into 



2 



Z^^{t) = Z{t) (A.4) 
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and, using the relation U^'^\t) = Z^'^\t)/ (l + Z^'^\t)\ Eq.([0|) becomes simply 



1 + Z(2)(t) 



(A.5) 



As in Section 4.2, in order to expand all the above generating functions in the 
scaling limit, it is convenient to treat separately their values right at ^ = 0, i.e. when 
9 — 9cmt — 1/12. Beyond the quantities Ai = ^^(1/12), A^-* = [Ai)'^ and their 
generating functions A{t) and A^'^\t) already defined in Sections 3.2 and 3.3, we also 
introduce the notations 



Af = J2 AjA-j , APP{t) = J2 = , A^^^PP{t) = J2 "^AAff . 

j=l i—l i—1 

(A.6) 

From Eqs.( |A.4D , p.lUD and ( |4.11D , we have the expansion 



from which we deduce the scaling behavior 



2nV4 
3s 



3 r 



(A.7) 



This expansion, together with the expansion (|4.9D for Zi — Ai leads to 

where we also made use of the large i leading values AijT' ~ 1/3 and 74^^/2* = 
n^/^r/Q. This yields 



^1/4 
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-1/4 



2nV4 
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(A.IO) 
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while, from 



(20 



3ni/4 
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(A.ll) 



we have 
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i(2) 
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+ ... (A.12) 



Inserting the expansions ( |A.1C1| ) and ( |A.12|) in (|A.5|) and using the leading behaviors 



^(2) j^g-sn ^/4/4j _ ni/V(9s) and i^^)?? (^e"^'^ '^V4) ~ 2nV2/(27s2), we obtain fi- 
naUy 
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n 



1/4 



3 r"^ 



+ . . 



upon integrations by part over r. Here we used the notation 



i(2)pp(i) _2 (i(2)(t) 



which is the critical value of V^'^^'^^if) at (7 = 1/12. We end up with 



+ ... 



(A.13) 



(A.14) 



1 _|^^(2)pp _ ^(2)pp^ 



(2^)2^^ 



n 



1/2 



24 



(A.15) 



where W^^"^^ is the coefficient of in W^'^2)pp(^^^ From the exact expression for A^, we 
can easily compute W^'^^'^'^if) exactly. In particular, we have the small r] expansion: 

yj7(2)pp _ = 4 + (12 _ 87r2 + 24 log(r7))r7 + 0{r]^) . (A.16) 

From the singularity at 77 = 0, we deduce that W^^"^^ jiT'Y ~ 24/z^ = 24/(7i^/^r^) in 
the scaling regime, and therefore 



1 



(2)pp 



(2^)2^ 

Upon integration over ^, this gives finally 



(A.17) 



V. 



(2)pp 



(2^) 



^^2 



n — >cx3 



IT 



Ttn 



3/2 



+00 



18 r / d^^e-^ 



(A.18) 
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Dividing this result by 



(2) 



in 



n — >co 



12' 



Txn 



7/4 



54 



+ 00 



yields the average number of contacts 



1/4 ^ 

3 



(A.19) 



(A.20) 



For small r, this gives the finite i result (c)^ = ^ of Eq.(3.26). The linear dependence 
of the average number of contacts extends in practice to the whole range of r in the 
scaling limit. 



Appendix B. Correction to the area correlation for two weakly avoiding 
geodesies 

The scaling result ( |4.28D , i.e. {ujiuj2)r = for two weakly avoiding confluent 
geodesies holds in the scaling limit ti ^ oo. We can easily get the flrst correction to scal- 
ing by computing the 1/n term in Eq.( [4.25| ) giving U^'^\g; Ai, A2)/(2*)^. In practice, as 
we eventually apply the operator l/(4(Ui/U2)9^/9/Ui9/U2, we need only the corresponding 
cross term depending on both and 112- This term reads 

81 J 92 

n 



^■(r, Hi)J^{r, 112) 



2 



J^{r - r',/i2) 



3r' J \ 3 (r — r') 



This leads immediately to the correction 

/ \ 9 1 r^^,, f^dj^ 

^ / r ^ , 1 a^. 1 



In particular, at small r, we have the expansion 



( r dr' — ^(r\a — ^(r-r\a 



(B.l) 



(B.2) 



i^^0 = 7^+^^^ + --- (B.3) 



2^ dC 30 ^140 

so that the { ■ } term in the integrand of Eq. (|B.2| ) behaves as 

100 ^350 ^ ^ 

Upon integration over ^, only the second term survives so that the integral behaves as 
i^/n X 9r^/700. From the small r behavior p(r) ~ 3r^/7, we deduce 



which is nothing but the announced result (|3.28|) . 
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